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Abstract. A mathematical model describing the initial stage of the capture into the 
parametric autoresonance in nonlinear oscillating systems with a dissipation is con¬ 
sidered. Solutions with unboundedly growing energy in time at infinity are associated 
with the autoresonance phenomenon. Stability of such solutions is investigated. We 
describe classes of admissible deterministic and random perturbations such that the 
stability of autoresonance is preserved on an asymptotically large interval. 
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1. Introduction 


Autoresonance is a phenomenon of a considerable growth of the energy of forced nonlinear 
systems. This phenomenon plays an important role in a wide class of various physical problems 
associated with nonlinear oscillations and waves [1-3]. The majority of well-known theoreti¬ 
cal investigations [4-9] consists a numerical and asymptotic analysis of mathematical models 
describing the initial stage of a capture into autoresonance. However, the stability of such 
models in the presence of external perturbations remained an open question. In the present 
work we study the problem of stability of parametric autoresonance in nonlinear systems under 
persistent perturbations. 

Consider a model system of primary parametric resonance equations [4]: 


( 1 ) 


dr . , r dxjj 

— — r sm ip — or, — 
dr dr 


r — Xr + fcos'i/j, r > 0. 


This system appears in the asymptotic analysis of nonlinear oscillations driven by a small force. 
The unknown functions r(r) and -0(r) represent the slow varying amplitude (energy) and phase 
shift of fast harmonic oscillations. The parameters A > 0 and f ^ 0 are factors related to the 
driving frequency and amplitude. The positive constant S correspond to dissipation coefficient. 
Solutions with unboundedly growing energy in time are associated with the capture of an 
oscillatory nonlinear system into parametric autoresonance. The aim of this paper to prove the 
stability of resonance solutions. It is assumed that only stable solutions correspond to motions 
that are observed in nature. Note that stability of the autoresonance in nondissipative systems 
(5 = 0) was discussed in [10]. 

System (1) is derived by an averaging of parametrically driven nonlinear oscillations [11]. 
Let us consider the equation 


(2) — + p— + (1 + £COS0) x + qx 3 = 0 

dt z dt 

as an example of the initial mathematical model. Here, 0(£; a) = 2 1 + at 2 , 0 < e, a, [3 -C 1, 
7 = const > 0. The point x = 0 is a stable equilibrium of the unperturbed oscillator (e = 0). 
Solutions of equation (2) with initial values near the equilibrium |x(0)| + |i(0)| <C 1 and with 
amplitudes increasing up to the order of unity at large times correspond to autoresonance (see 
Figure 1). 
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FIGURE 1. The energy of the parametrically driven nonlinear oscillator (2). Pa¬ 
rameter values are /? = 0, 7 = 3/2, £ = 0.001, a = £ 2 /8. Two curves correspond 
to different initial points. 


The asymptotic of resonance solutions of equation (2) for 0 < t <C £ 2 is constructed in the 
form: 

^ £t 2 

x(t ) = a/ k£t{t) cos -(0 + 0(r)) + 0 (e), £ -► 0, r = —, « = —, 

where the slow varying functions r(r) and 0(r) satisfy system (1) with A = 8cc£ -2 , 5 = 2(3e~ 1 , 
and f — 1. 


2. Autoresonance solutions 

The solutions of system (1) cannot be written in an explicit form. However, the asymptotic 
solutions with increased energy at infinity r —> 00 can be constructed in the form of power 
asymptotic series with constant coefficients: 


( 3 ) 


R±(r) = Xt + Y^ 
3=0 


r± t ^ 
3 ’ 




T = 


3 =0 


± T~ j 
3 1 ’ 


t —y 00 . 


Substituting these series in system (1) and equating the expressions of the same powers give 
the recurrence relations for determining the coefficients rf and In this way two solutions 


are constructed if 0 
sin = S: 


< 5 < 1; difference stem from two roots of the trigonometric equation 


Tpf = 


COS -0, 


± > 


rt = ~f cos0^, rf = /tan0^. 


± 




Series (3) correspond to the exact solutions of system (1) with a given asymptotic expansion at 
infinity [12]. We investigate the stability of the solutions R±(t), \I/±(t) in the sense of Lyapunov 
and under persistent perturbations. 

The solution with the phase ^ = arcsin S is unstable as can be seen by analyzing the 
equations linearized near the leading term of the asymptotic solution: one of the characteristic 
roots has a positive real part. In the case of 00 = n — arcsin S, such an approach is unapplicable 
because the characteristic roots are purely imaginary. In this situation the property of stability 
depends on non-linear and time-dependent terms of equations (see [13,14]). 


3. Perturbed equations 
Along with (1), we consider the perturbed system 


( 4 ) 


dv d'lp 

— = (1 + fi£)r sin 0 — Sr, — = r — At + pC + (/ + nrj) cos 0, 
dr dr 
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where the perturbations £(r, ip,r), p(r,-0 ,t), and ((r,if>,T) are defined for G R 2 , r > 0. 

The coefficient /t 6 1, 0 < /i C 1 is a perturbation parameter. We consider the perturbations 
such that system (4) has a global solution. To ensure this property we have to require additional 
restrictions (see [15,16]) on the class of functions (£,p, C)- Our goal is to identify a class of 
persistent perturbations V such that any solution t m (t), Vv( r ) °f system (4) with initial data 
from a neighborhood of (7) remains near the solution T_(t) while the parameter p is 

small and (£, 77 , C) belong to V. 

The example of original system leading to (4) is the Duffing oscillator 

1 + e(l + /m) cos(</> + 7793 ) | x + yx 3 = 0 , 0 < /i < 1 , 

where the functions a(x,x,t]e) and tp(x,x,t]e) correspond to perturbations of the pumping 
amplitude and phase. If a = aft ; e) and <p = <p(t; e), then the perturbations in system (4) have 
the form 

(6) £(T) = a(t;s), tj(t) = a(t;e), ((r) = - 4cPt ^ t,£ \ t =—. 



4. Lyapunov stability 


To study the stability of the solution R_(r), 4te(r), we use the first terms of asymptotic 
expansion (3) 

(7) i?_(r) = At +/cr + C>(r _1 ), T_(r) = 7r — arcsin 5 — — t _1 + 


a 


where a = \/l — 5 2 . 


Theorem 1. If 0 < S < 1 and f > 0, then the solution 4te(r) with asymptotics (7) is 

asymptotically stable. 

Proof. By the change of variables 

(8) r = R-(r) + VXr R, V’ = ^ 

system (1) can be rewritten in the form 
1 


Vx t dT 


dR = -BtH(R,<S!,T) + F(R,<S!,T), -L If=d R H(R,V,T), 


( 9 ) 

where 

Jj>2 

H(R, T, r) = — + — 1 cos(T + T_) — cos 4te 
2 At . 


VXr dr 


F(R, T, t) = 


R 


\/At - 


5 + (/ — 1) sin(T + 4te) 


T sin 
R 


fR 

yf \r 


cos(T + T_) — cos 


2 T 


-1 


For the new functions R{r ), T(t) we study the problem of stability of the equilibrium (0; 0) by 
the Lyapunov second method. To construct a Lyapunov function for system (9) the asymptotics 
of the right-hand sides in a neighborhood of the equilibrium (as p = yjR 2 + T 2 — > 0) and at 
infinity (as t —> 00 ) are used. Note that all asymptotic estimates written out bellow in the 
form 0(p n ) and 0(T~ m ) (n, m = const > 0) are uniform with respect to R, T, t in the domain 

B(p *, t*) = {(i?, T, t) : p < p*, t > t*}, p*, t* = const > 0. 

It can easily be checked that the Hamiltonian has a positive quadratic form as the leading 
term of the asymptotic expansion: 

H = ^- + a^- + 0{p 3 ) + 0{p 2 )0(T~ 1/2 ), p^0, t y 00 . 
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By taking into account (7) one can readily write out asymptotics of the derivatives: 

d R H = R-\ — j= [cr(l — cosT) — 5 sin 'hi + 0 (p)C>(t“ 3 / 2 ), 

V At 

d^H = a sin \l/ + <5(1 — cos \P)4— ^L= \a sin T + <5(1 — cos T)1 + 0(p)0(r~ 1 ), 

V At 


d r H = 0( P 2 )0(t~V 2 ). 

The non-Hamiltonian part F(R, T, t) tends to zero as t —> oo: 

F =-m[R + 0(p 2 )]r- 1/2 + 0( P )0 (t- 1 ), m 


K 

y/X 


> 0. 


The Lyapunov function is constructed on the basis of the Hamiltonian: 

(10) V(R,V,t) = H(R,^,t) + 7 ^-R^t- 1/2 . 

The derivative of the function V(R,^,t) along the trajectories of system (9) decreases as 
t —y oo; the leading term of its asymptotic expansion consists a quadratic form: 


1 dV 


Vx t dr 


(9) 


1 dV dV\_dH_ - 

i/At dr dR . . 


dV dH 

OR 


= [R 2 + a T 2 ] [1 + O(p) + Oir- 1 )}^ 1 / 2 . 

Note that the remainders can be made arbitrarily small by choosing suitable domain H(p*,T*). 
It follows that there exist pi > 0 and T\ > 0 such that inequality 


dV 

dr 


(9) 


?77/\'A / O ox 

<- (R 2 + aT 2 ) 


holds for any (R, T,t) G £>(pi,Ti). Similarly, there exists p 2 > 0, t 2 > 0 such that 

(11) \{R 2 + aT 2 ) < V(R, T, t) < \{R 2 + aT 2 ), V (R, T, t) G B(p 2l r 2 ). 

Thus 


( 12 ) 


dV < my/X 
dr (9) ~ 3 


for each triple (R, \P, r) G B(po,To), where po = min{pi,p 2 } and t 0 = max{Ti,T 2 }. Let e be an 
arbitrary positive constant such that 0 < e < po; then 

(13) sup V(R, T, t) < — 1 - < —— < inf V(R, \P,r), 5 e = e 

p<S e ,T>r 0 4 4 p=€,r>r 0 



Therefore any solution R(t), T(t) of system (9) with initial data [-R 2 (to)+'L 2 (t 0 )] 1 / 2 < 5 e cannot 
leave e-neighborhood of the equilibrium (0; 0) as t > t 0 : [R 2 (t) + T 2 (t)] 1//2 < e. Integrating 

(12) with respect to t, we obtain 


0 <V(R, T, t) < C exp(—2 /t) V(E,Lt) G%,r 0 ), 

where positive constant C depends on a trajectory of system (9), l = m\/X/§ > 0. Hence 
the Lyapunov function tends exponentially to zero along the trajectories of system (9). If we 
combine this with (8) and (11), we obtain asymptotic estimates for solutions of system (1) with 
initial data from a neighborhood of the solution R_(r ), T_(t): 

r(r) = R-{t) + 0(T 1/2 e~ lT ), ^{t) = T_(t) + 0(e~ lT ), t > t 0 . 
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This completes the proof. 

Actually, the sufficient conditions obtained in Theorem 1 are almost necessary; this can be 
seen from the following theorem. 

Theorem 2. If 0 < 5 < 1 and f <0, then the solution 7?_(r), \I/_(t) with asymptotics (7) is 
unstable. 


The proof is based on Lyapunov function (10). Since / < 0, we have 


dV_ 

dr 


> + a^ 2 ) > 0 

(9) - 4 


for any (R , T, r) G £>(po, fi))- From the Lyapunov theorem (see [17]) it follows that the solution 
i?„(r), T_(t) is unstable. 


5. Deterministic perturbations 

In this section we consider the problem of stability of the capture into the parametric autores¬ 
onance under persistent deterministic perturbations. Our goal is to identify a class of functions 
such that system (4) has resonance solutions with growing energy as the perturbation 
parameter p is sufficiently small. 

Let 7/ be a positive function such that 7/ —> oo as p —> 0. Now we shall give the following 
definition of stability (see [18]). 

Definition 1. The solution R_(r), T_(r) of system (1) is stable under persistent perturbations 
V on an asymptotically large interval (r 0 ; r 0 + 7/) if Ve > 0 3 <5 e , A e > 0 : 

V Qo, 0o : |0o - #-( T o)| + |0o - 'M r o)| < <5e, Vp< A e , V(£, r?, () G 7^ 

the solution r^{r), ^(r) of system (4) with initial data r M (T 0 ) = Qo, Vv( T o) = 4>o satisfies the 
inequality 

sup |r M (r) - A_(t)|t -1/2 + |0 m (t) - ^_(r)| < e. 

0<r—ro<7^ 

Note that this definition is different from a classical one [17] because of the finite time interval. 
But such an approach seems to be reasonable, since the considered mathematical model (1) is 
valid only for 0 < r e _1 , where 0 < e <C 1 is the driving amplitude in the nonlinear systems 
like (2). In order to describe the autoresonance phenomenon for r e^ 1 one should consider 
other equations [5]. 

Consider a class V a ^ c of functions (£, 77 , </) such that £ T^a,b,c 

sup \f(r,'ijj,T)\T~ a + |?7(r,V’,T)|T -6 + |C(f,- 0, r)|r _c < oo. 

(r,i/>)GR 2 ,T>0 

Let h > 0 be a positive constant; we define a class I);;', r as a subset of such that for any 

e 

sup |£(r,0,T)|T“ a + |r/(r,'0,T)|T- b + |C(r, 0, t)|t“ c < h. 

(r,'!/))GR 2 ,r>0 

The Cauchy problem for perturbed system (4) with initial data from a neighborhood of the 
solution 72_(r), T_(t) is assumed to have a global solution. This requirement impose the 
additional restrictions on the class of perturbations (see [15]). 

Theorem 3. If 0 < 6 < 1, / > 0, then Vh > 0, a > —1/2, b > 0, c > 0, and x G (0; x 0 ) 
the solution R-(r), T_(r) with asymptotics (7) is stable under the persistent perturbations 

K.tf.C) £ Kac- 

sup |5(r,'0,T)|r - “ + \T)(r ) tp,r)\r~ b + |C(<W,t)|t-- c < h 

(r,V')GR 2 ,r>0 
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on an asymptotically large interval 0 < r < 0 (p x ), where x 0 = id 1 , id = max{a + 1/2, b, c}. 


Proof. As above, we reduce the considered problem to the analysis of the equilibrium (0; 0) 
in system (9). By change of variables (8) perturbed system (4) is reduced to the differential 
equations 

(14) -yL ^ = -d*H + F + pG, -^=^- = d R H + pQ, 

VAr dr VAr dr 

where 

= (R_ + RV *r)sm(tt + tf_)-^-, Q0R,T,r) = (/t)cos(^ + tf_) + C)^L. 

At V At 

Persistent perturbations of system (9) are associated with the functions G and Q. The functions 
£(R, T,t), rj(R, dqr), and </(/?, T,t) are associated with £(r, 0, r), p(r,ijj,r), and ((r,ijj,r) 
through (8), e.g., £(!?, T, t) = £(1?_(t) + Ry/Xr, T_(t) + T, r). 

Let h > 0, a > —1/2, b > 0, and c > 0 be arbitrary constants. The derivative of the 
Lyapunov function (10) with respect to t along the trajectories of perturbed system (14) has 
the form: 


(15) 


dV_ 

dr 


(14) 


dV 

dr 


+ pV\ r{Gd R V + Qdn,V). 

(9) 


Note that the first term in the right-hand side of (15) satisfies the inequality (12) in the domain 
B(p 0 , t 0 ), while the derivatives d R V, d^V are bounded: \d R V\ + \d^V\ < i. From the definition 
of the class T>% b c it follows that 


vT((G| + (Cl) < M h T*, 6 B(po,t 0 ), 


where Mh is the positive constant and id > 0. Thus the inequality 


dV 

dr 


< 


m 


(14) 


\/Ar 6 £M h £ 
V — n - r 


m 


holds for any ( R , T, r) G B(po, To). For any e and x such that 0 < e < po and 0 < x < l/id we 
define 




^ = <A/-, A e = 
o 


crh/m 


2A(2r o ym h 


l/z 


Then the derivative of the function V is negative: 

dV_ 
dr 


^ mV A 

"N 

1 

1 ^ 

On 

bO 

(14) “ 3 

L 4 J 


z = 1 — idx > 0. 


< 0 


if p < A e , S e < p < po and 0 < r — To < rop x . Let us remember that the Lyapunov 
function satisfies inequalities (13). Hence any solution R^r), v h A1 (r) of perturbed system (14) 
such that [1? 2 (to) + d* 2 (to)] 1//2 < 5 e cannot leave the e-neighborhood of the equilibrium (0; 0) as 
0 < t — t 0 < r 0 p~ x . Taking into account (8), we obtain 

Mt) - R— (t)|(At) -1/2 + |-0 m (t) - T_(r)| < e, 0 < r - r 0 < 0{p~ x ). 

From [15] it follows that the solution R_(r), 'L^(r) is stable on the finite interval (0; Tq]. 
Therefore for any h > 0, a > —1/2, b > 0, c > 0, and x G (0; x 0 ) the solution f2_(r), T_(t) 
of system (1) is stable under the persistent perturbations (£, p, () G V h a h c on an asymptotically 
large interval 0 < t < (P(/i - ^). 


Theorem 4. If 0 < 5 < 1, f > 0, then V h > 0, a < —1/2, b < 0, and c < 0 the solution 
R_(r), ^_(t) with asymptotics (7) is stable under persistent perturbations (£,p,C) £ : 

sup |£(r,'0,T)|T“ a + |r/(r,^,T)|T -6 + |C(r, -0, t)|t" c < h 
(r,V')eK 2 ,r>0 
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on the infinite interval r > 0. 

The proof is follows from Malkin’s theorem [17]. 

6 . Random perturbations 
W e consider the perturbed system: 

(16) — = (1 + f)r sin <f> — dr, — = r — At + C + (/+ 77 ) cost/;, 

dr ' dr 

where £(r, if, r; u, p), rj(r,i/;,T]UJ, p), ((r, r; uj, /i) are one-dimensional random processes de¬ 
fined on a probability space (Q, T. P), (r, ff) e M 2 , r > 0, p > 0. Our goal is to find a class of 
perturbations (£, 77 , () such that the stability of autoresonance is preserved in perturbed system 
(16). 

Let % be a positive function such that 7) t —> oo as p —> 0. Now we give the following 
definition of stability under random perturbations. 

Definition 2. The solution i?_(r), 4'_(r) of system (1) is stable under random perturbations 
V on an asymptotically large interval (r 0 ; r 0 + Tfi) if Ve, v > 0 3 S £ , A > 0 : 

V go, 0o : 1 00 - 7?_(t 0 )| + |-0o - df_(r 0 )| < A, V/x < A, V(^, 77 , ()6P 

t/ie solution r^fir^oj), of system (16) with initial data r^T^uj) = g 0 , = 0 O 

satisfies the inequality 

p( sup jr A 1 (r;o;) - R_(t)|t _1/2 + |^(r;a;) - T_(r)| > e) < n. 

' 0<r—ro<7^i ' 

Note that this definition is usually used for strong stability (see [16, p. 152], [18, Chap. 9], 
and [19, p. 400]). 

Consider a class 7 Z a ^ c of random functions £, 77 , (. Assume that for all (£, 77 , £) 6 'Tt a .b, c there 
exists at least one random function S(r-,uj,p) such that 

T+l 

3 1 /( 0 ;) > 0 : M r S d = j S(t',u, p) dt < pv(uj) Vr> 0, uj E Cl, p > 0, 

T 

and 

Ez/ l = j tt(u;)P (dcu) < 00 . 
n 

It is assumed that for all (£, 77 , £) G 77 ai ;, )C 

sup |£(r-,' 0 , 7 -;a;, 77 )|T _a + ^(r, ip, r; cu, /i)|r -6 + |C(r,-0, r; tu,/r)|T“ c < ^(rja;) 

(r,^)gR 2 

for all r > 0 , a; G fl, and p > 0 . 

Let h > 0 be a positive constant; we define a class 77(] 6 c as a subset of 77 aj b iC such that 
E v < h. 

Notice that the stability under random perturbations with bounded expectation (sup t El? < 
00 ) was investigated in [16, p. 26] provided the unperturbed system is dissipative in the sense 
of [16, p. 8 ]. However, the considered equations do not have this property because system (1) 
has solutions of two types: with bounded and unlimited amplitudes, see Figure 2. Thus the 
results of [16] are not applicable. Stability of nondissipative systems is discussed in [10,20]. We 
extend these results to the analysis of stability on an asymptotically large interval. 
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Figure 2. System (1) has solutions with bounded and unlimited amplitude. 
Two curves correspond to different initial data: (left) r(0) = 1.59, -0(0) = 0.59, 
(right) r(0) = 0.35, 0(0) = 3.09. Parameter values are / = 0 .2, S = 0.5, and 
A = 1. 


Theorem 5. If 0 < 5 < 1, / > 0, then V h > 0, a > —1/2, b > 0, c > 0, and x 6 (0; x 0 ) 
the solution 7?_(t ), \I/_(t) with asymptotics (7) is stable under random perturbations (£,??, C) G 
V h ■ 

'*'a,b,c* 

sup \^\r~ a + \p\r^ b + \(\t~ c < S(t;uj, p) Vr> 0, u; G ff, p>0 

(r,ip)£ K 2 

sup Mt-S 1 < p z/(cu), E// < h 

T >0 

on an asymptotically large interval 0 < r < 0(p~ x ), where x 0 = d~ l , d = max{a + 1/2, b, c}. 


Proof. We reduce the problem to the analysis of the equilibrium (0; 0) in system (9). Change 
of variables (8) leads to the perturbed system 

(17) -^=^ = -dyH + F + G, -}=^- = d R H + Q, 

V At dr v At dr 

where random functions G(R , T, t; uj , p), Q(R , T, t; uj , p) are defined as follows: 


G = (R- + rV\t) sin(T + T_)-^, Q = (/? 7 cos(T + T_) + £) — 1=- 

At vA t 

Let e > 0, v > 0, h > 0, a > — 1/2, b > 0, and c > 0 be arbitrary constants. We construct 
a Lyapunov function for system (17) on the basis of the Lyapunov function V(R, T, t) for the 
unperturbed system (see [21]): 

U{R, T, t; uj, p) = V(R, T, t) exp $(t; uj, p), 

where a smooth function $(t;u;,/ x) is defined bellow. The derivative of U(R,^,t]uj, p) with 
respect to t along the trajectories of perturbed system (17) has the form: 


dU 

dr 


= U + (+ V\r(G d R V + Q cApF)) exp d>. 
V dr ( 9 ) / 


(17) 


From the definition of the class 77/ h c it follows that 


(|G| + IQI)r 1/2 -'’ (fl,*,r) 6 B(p„,r„), 


where q is the positive constant and i7 > 0. The partial derivatives d R V , c^F satisfy the 
inequality: |<9 rF| + \d^V\ < p£/po in the domain B(p 0 ,r 0 ). Hence in view of (11) and (12) the 
derivative of U(R, T, t; uj , p) satisfy the estimate: 


dU 

dr 


(17) 


< d^U - 


m 




1-pR 


24 £q S 
ma5 e p 0 p 


U 
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in the annular domain 5 e < p < e, r > r 0 . For any x such that 0 < x< 1/d we define a 
random variable 


A„, = 


Ap 

v(uj) 


An — 


mcr5 e po 
24 (2r 0 )*V 


z = 


1 — dx 


> 0. 


Then for any /i < A w the inequality 

dU_ 
dr 


< driU _^( v - s _' w 
(17) 3 v V 




holds in the domain < p < e, 0 < r — r 0 < r 0 /i x . 
For a fixed ca G 12 we consider the integral 

k+1 

S(t\u,p) 


I(k;u,fi) = / u(u) - 


p 


dt, k = 0,1,2,.... 


From the definition of the class 7 Z l / bc it follows that I(k] 0 J,p) > 0 for any k > 0 and ca G 12. 
Define an auxiliary random function 9(t]0J,p) such that 

k+l 


(18) 


9(t;cu,p)dt = I(k;oj,p) V k > 0. 


Since right-hand side of (18) is not negative, it follows that there exists a non-negative function 
9(t]u, p) > 0. Without loss of generality, we can assume that 0(r; oj, p) is a continuous function 
such that 9(k] ui, p) — 0 for any k = 0,1,2,.... Let us define <f>(r; co) as follows 




m 


Vx 


/ \ ta, /i) { ^ 

z^(cc)-d(t; ca, p) dt. 


p 


Then the derivative of U satisfies the inequality: 

dU 
dr 


<-HRR,u< o 

(17) 3 V 


in the domain 5 e < p < e, 0 < r — r 0 < r 0 p _x . Taking into account properties of the functions 
9(t]u,p) and S(t;oj,p), we obtain |<L(r;a;,/i)| < $ 0 , $0 = 4m\/A/3. Thus for any e > 0, 
u) G 12, and p > 0 we have 

sup !/(/?, T, r; (u,/i) < ^-exp $0 < exp(—<f>o) < inf U(R,^,r;u, p), 

p<S e ,T>T 0 4 4 P=£,T>T 0 

where S e = e exp(—<&o)y/cr/6. Consider perturbations (<C, 7 7, 0 £ hZabc suc h that u ( u ) — h/ v 
uniformly for all uj G 12. Then the parameter p can be bounded away from zero 0 < p < A = 
(i>A o/h) z < A u . Therefore any solution i?(r;u;), T(r;cj) starting from the neighborhood of 


equilibrium [i? 2 (r 0 ; u) + T 2 (r 0 ; <u)] 1/2 < S e 


remains 


inside the ball [f? 2 (r; u) + T 2 (r; u;)] 1 / 2 < 


for 0 < t — tq < 0(p~ K ). 

For all other perturbations such that v{uj) > h/v, it follows from the Chebyshev inequality 
that 


P(z/(ca) > h/v) < 


E v 

h/v 


< v. 


In this case, solutions of perturbed system (14) can leave any e-neighbor hood of the equilibrium 
(0; 0). However, the probability of such events is small: 


p ( sup [ R l( r 'iu) + ^ 2 (t;uj )} 1/2 > e) < v, % = 


0<r—ro<7)i 


T 0 /l 
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The change-of-variables formula (8) implies the stability of the solution i?_(r), T_(r) of 
system (1) under random perturbations on the asymptotically large interval 0 < r — To < 
0{p~ x ). Stability on the interval r G (0; To] follows from the continuity of solutions with 
respect to parameters, see [16]. This concludes the proof. 

Theorem 6. If 0 < 5 < 1, / > 0, then V h > 0, a < —1/2, b < 0, and c < 0 the solution 
R-(t), \I/_(t) with asymptotics (7) is stable under random perturbations (£, 77 , C) £ 77 % bc '■ 

sup \^\r~ a + \r]\r~ b + \(\r~ c < S(t;cu, fi) Vr> 0, to G fi, /1 > 0, 

(r,'0)eR 2 

sup M r S' < /iz/( u>), E v < h 

T> 0 

on the infinite interval r > 0. 

The proof follows from the results of [20] and [10]. 


7. Examples 


1. To illustrate Theorem 3, let us consider the following example. Let 0 < S < 1, / > 0, 
£(r) = 1, r](r) = 1, and £( t ) = r. Then perturbed system (4) takes the form: 


dr 

dr 


(1 + /i)r sin 7 /; — 5r, 



At + /rr + (/ + / 1 ) cos "0, 


0 < /i <C 1. 


Note that one can easily construct the asymptotic expansion at infinity r —> 00 for a particular 
resonant solution of the perturbed system as follows 


c 

r u( T ) = (^ - h) T + 0 ( 1 ), Vv( r ) = n - arcsin (yyp^) + 0 (r _1 ). 
Therefor we have 


|7?-(t) — t m (t)|t 1/2 = /ir 1/2 [l + 0(t 1 )], |^_(t) - ^(t)| = O(p) 

as t —y 00 and /i — > 0. From the first estimate it follows that the solution AL(t), T_(t) is 
unstable for r > /i -2 . However, since £, 77 , £ belong to 27 >q 015 we see that from Theorem 3 
it follows that the solution R_(t), T_(t) is stable under the persistent perturbations at the 
interval 0 < r < 0 (p~ x ), where 0 < x < 1 . 

2. Let us consider random perturbations of the prime parametric resonance equations. First 
define the random process 

N 

J n {t;uj,h) = £ j n (co)x(n < t < n + /i), 0 < /i <g; 1, r > 0, uGO, JVgN, 

72—1 

where y is a characteristic function, {j n (co)}^ =1 are random variables such that E \j n \ < 00 . 
Since M r |J/v| < p[\j n \ + |j n +i|] for T £ [n; n + 1), we see that the random functions 

£(r; u, p) = w, h), r/(-r; cu, /1) = Jtv(t; w, 77), £(t; u, fi) = tJjv(t; w, //), 

belong to the class 77. 0 i o,i- From theorem 5 it follows that the solution R_(r), T_(r) is stable 
as r < /i _1 . 

3. Another example of perturbations is described by the function with a random jump. Let 

£(t;w,/x) = J(r;u;,/i), 7?(r; cu,/x) = tJ(t; cu,/r), (i), 

where J(t; cu, / 1 ) = j(uj)x(w < t < a; + / 1 ), 0 < /i 1, r > 0, u G fl = (0, 00 ), and j(iv) ^ 0 is 
a random variable with the bounded expectation E|j| < 00 . Since M r | J(r;cu,/i)| < /i|j(a;)| for 
r > 0, it follows that (£, 77 , £) G 77. 0 i i,i and the solution AL(t), T_(t) is stable as r <C /i” 1 . 
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Finally note that from Theorem 3 and (6) it follows that stability of the capture into au¬ 
toresonance in perturbed system (5) is preserved on the interval 0 < et <C p~ l if 

0 </3<^, \a(t;e)\ < fiy/et, \<p(t; e)| < pE 2 t 2 , t > 0. 

8. Conclusion 

Stability of the capture into the parametric autoresonance in dissipative systems under per¬ 
sistent perturbations on a long time interval was proved. The length of the interval depends 
on both the perturbation parameter /i and the class V of admissible perturbations. Stability 
of autoresonance under white noise perturbations was not considered. This will be discussed 
in a further paper. 
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